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Abstract 

Based on empirical market data, a stochastic volatility model is 
proposed with volatility driven by fractional noise. The model is used 
to obtain a risk-neutrality option pricing formula and an option pricing 
equation. 
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1 Introduction 

The Black-Scholes (B-S) £Q [2] pricing formula for call options is widely used 
by market agents. However, rather than as a predictive formula, it is used as 
a means of quoting option prices in terms of another parameter, the implied 
volatility. If Cbs (St, K,T — t, a t ) is the Black-Scholes price for an option at 
time t, with strike price K, maturity time T and underlying price St then, 
the implied volatility <j l t mp is such that 

C BS (S t ,K,T-t,aD =C t (K,T) (1) 

C t (K, T) being the market price of a (liquid) option. 
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B-S is based on unrealistic assumptions about the market process: geo- 
metric Brownian motion with constant non-stochastic volatility, continuous 
adjustment of the portfolio and no transaction fees. Nevertheless B-S owes 
its great popularity to its simplicity and to the fact that Cbs is an invertible 
map from [0, oo) to [0, St — Ke~ r ^ T ~ t ^ , r being the risk-free rate. 

The empirical shortcomings of B-S originated a vast amount of alterna- 
tive models j3] — which attempt to explain the deviations from B-S by 
introducing additional degrees of freedom. The majority of these models 
have no compact closed-form solution and numerical solutions most often do 
not reproduce correctly the data profiles (TU] [TT] . 

Given the success of the Black-Scholes formula as a compact invertible 
parametrization, any alternative model that attempts to introduce a higher 
degree of market realism should, at least, provide also a compact closed-form 
recipe to parametrize the data. 

In liquid markets the autocorrelation of price changes decays to negli- 
gible quantities in a few minutes, consistent with the absence of long term 
statistical arbitrage. Geometric Brownian motion models well this lack of 
memory, although it does not reproduce the empirical leptokurtosis. On the 
other hand, nonlinear functions of the returns exhibit significant positive au- 
tocorrelation. For example, there is volatility clustering, with large returns 
expected to be followed by large returns and small returns by small returns 
(of either sign). This has the clear implication that, on stochastic volatility 
models, long memory effects should be represented in the volatility process. 

Two-factor models with one persistent factor and one quickly mean- 
reverting factor ^2] > multicomponent GARCH models as the one proposed 
in (LM(q)-ARCH) or fractionally integrated processes [H] [T^] provide 
adequate fits to the empirical volatility data. However, these models do 
not possess the analytical properties required to construct an option pricing 
equation generalizing Black-Scholes. 

Because the market process, as well as some functionals of their variables, 
display approximate self-similar properties, mathematical simplicity suggests 
to look for descriptions in terms of fractional Brownian motion. In fact, if 
a nondegenerate process X t has finite variance, stationary increments and is 
self-similar 

Law (X at ) = Law (a 11 X t ) (2) 
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then jlH] < H < 1 and 

Cov (X s , X t ) = \{\s\ 2H + \t\ 2H -\s- t\ 2H ) E {Xl) (3) 

The simplest process with these properties is a Gaussian process called frac- 
tional Brownian motion. Fractional Brownian motion J7| 

E [B H it)] =0 E [B H (t) B H (s)} = l - { \t\ 2H + \s\ 2H -\t- s\ 2H ] (4) 

has for H > | a long range dependence 

oo 

Cov (B H (1) ,B H (n + l)- B H (n)) = oo (5) 

n=l 

and was suggested in the past as a tool for modeling in Finance (THJ - However, 
because it was pointed out [TH] that markets based on Bh it) could have 
arbitrage, fractional Brownian motion was no longer considered, by many, as 
promising for mathematical modeling in Finance. The arbitrage result in 19 
is a consequence of using pathwise integration. With a different definition 

en, 

/ / (t, u) dB H (t) = lim V / (i fc , w) (B H (*fc+i) - B H (t k )) (6) 

where A : a = to < ^i < • ■ ■ < t n = b is a partition of the interval [a, b], 
|A| = maxo<K„_i (4+1 — tfc) and o denotes the Wick product, the integral 
has zero expectation value and the arbitrage result is no longer true. This 
is, in fact, the most natural definition because it is the Wick product that is 
associated to integrals of Ito type, whereas the usual product is natural for 
integrals of Stratonovich type. 

An essentially equivalent approach constructs the stochastic integral through 
the divergence operator and Malliavin calculus 22 . A fully consistent stochas- 
tic calculus has therefore been developed for fractional Brownian motion [211] 
fH\ j2Sj [23] [22] |2f>j . We will draw on these results, in particular on the frac- 
tional generalization of the Ito formula, to derive our stochastic differential 
market model and option pricing equation. 

However, to simply postulate that the price process follows a stochastic 

law 

dS t = fiS t dt + aS t dB H (t) (7) 
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with Bh (t) (H > 1/2) a persistent fractional Brownian motion, as done in 
EHj EZ], although mathematically pleasant, is not phenomeno logical correct 
in view of the fact that price changes have a short memory. The simplest 
alternative would be to generalize current stochastic volatility models by 
changing the Hurst coefficient of the volatility process, namely 

dS t = fiStdt + cr t S t dB (t) , . 

da t = iiy (9 — cr t ) dt + GyOtdBn (t) 

The price would follow a geometrical Brownian process, whereas the volatility 
would have a mean-reverting deterministic term and a stochastic component 
modeled by a geometrical fractional Brownian motion of Hurst coefficient 
H > |. This, on the one hand, would introduce a memory component on 
the volatility (where the data says it should be) and, on the other hand, 
as it is known [2H] (HHl stochastic volatility changes the effective probability 
distribution function of St (even for H = |), generating fat tails closer to the 
empirical data. 

However, the stochastic volatility model of Eqs.(jSJ), as we will see later, 
turns out not to be justified as well, because its scaling properties are very 
different from those of the empirical data. For this reason we have made a 
detailed analysis of the data, reported in Section 2, trying to infer a model 
which is, at the same time, mathematically manageable and not grossly dis- 
proved by the data. Our conclusion is the following coupled stochastic system 

dS t = fiS t dt + a t S t dB (t) ( 
\oga t = p+k{B H (t)-B H (t-5)} W 

It means that, in addition to a mean value, volatility is driven not by frac- 
tional Brownian motion but by fractional noise. Notice that our empirically 
based model is quite different from the usual stochastic volatility models 
which assume the volatility to follow an arithmetic or geometric Brownian 
process. 5 is the observation scale of the process. In the 5 — > limit the 
driving process would be the distribution-valued process Wh 

W h = lim - 5 (B H (t) — Bjj (t — 5)) (10) 

In Q the constant k measures the strength of the volatility randomness. 
Although phenomenologically grounded and mathematically well specified, 
the stochastic system (j^J) is still a limited model because, in particular, the 
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fact that the volatility is not correlated with the price process excludes the 
modeling of leverage effects. 

In Section 3, following a risk- neutrality approach jl] [Zj, we use the con- 
ditional probability distributions following from Eqs.tjHJ) to derive an option 
price formula and compare it with Black-Scholes, namely in terms of the 
equivalent implied volatility surfaces. 

The risk-neutrality approach of Section 3 provides a fairly simple explicit 
generalization of Black-Scholes, which we believe to be based in a more realis- 
tic mathematical representation of the market process. Nevertheless it relies 
on some approximations which may or may not be justified. Therefore, for 
future reference, we derive in Section 4 an option pricing equation of wider 
generality and obtain an integral representation of its solution. 



2 The volatility process. An empirical anal- 
ysis 

Option pricing being our primary motivation, the relevant time scales and 
temporal horizon is greater than one day. Therefore we will use daily data 
from the New York Stock Exchange (NYSE), the aggregate index and indi- 
vidual companies as well. To discount trend effects and approach asymptotic 
stationarity of the processes, we have detrended and rescaled the data as ex- 
plained in Ref. |31j. 

The first objective is to check whether the assumption that the price 
process follows a geometric Brownian process 

dS t = S t (fidt + adB t ) (11) 

is an acceptable hypothesis for the daily data. For detrended data \i = and 
if the stochastic part is a self-similar process it should satisfy 



E 



S (t + A) - S (t) 



S(t) 



A H (12) 



In the lower left panel of Fig.l one checks Eq.(|12p for the NYSE index 
in the period 1966—2000 with A = 1 : 30 days. We obtain approximate 
linearity in the log-log plot with H = 0.52. There is a small deviation from 
linearity for the first few days, but one should notice that stochastic volatility 
effects have not yet been taken into account. So, the geometric Brownian 
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Figure 1: Self-similar properties of the price S (t) and integrated log- volatility 

(3t + R a (t) processes 



hypothesis for the daily price process is, at least, not very unreasonable. It 
would imply that, in first approximation, the probability distribution would 
be log-normal with volatility (variance) given by 



\T - Ti 



■var (log S t ) 



(13) 



Eq.ffTB^) is used to extract the volatility process a t from the data and one 
checks whether a relation of the form 



or 



E\a(t + A) -a(t)\ ~ A H 
a(t + A)-a (t) 



E 



a(t) 



A 



H 



(14) 



(15) 



hold for the volatility process. This would be the behavior implied by most 
of the stochastic volatility models that have been proposed in the past. It 
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turns out that the data shows this to be a very bad hypothesis, meaning 
that the volatility process itself is not self-similar. However, the integrated 
log-volatility is well represented by a relation of the form 

t/s 

J2^gcr(n5) =(3t + R a (t) (16) 

n=0 

As shown in the lower right panel of Fig.l the R a (t) process has self-similar 
properties. We will identify it with a fractional Brownian motion with Hurst 
coefficient H ~ 0.8 (for the NYSE index), 

R a (t) = kB H (t) (17) 

The same parametrization holds for the data of all individual companies 
that we tested, with H in the range 0.8 — 0.9. This analysis bears some 
resemblance to the detrended fluctuation analysis technique [22]. Notice 
however that here, because we are using market data that has already been 
detrended and rescaled, it is used merely to extract the mean volatility. 

Eq. (jl6|) provides a simple mathematical representation of the volatility 
memory effects. In particular, it means that the volatility is not driven by 
fractional Brownian motion but by fractional noise, 

log<r t = (3 + ~ (B H (t) -B H (t- 5)) (18) 



S being the observation time scale (one day, for daily data). This provides 
a simple interpretation of the fact that empirical return statistics depend on 
the observation time scale. In the 5 — > limit (continuous time resolution) 
the driving process would be the distribution valued process (fractional white 
noise) 

W H (t) = lim - 5 (B H (t) — Bjj (t — 5)) (19) 
For the volatility (at resolution S) 

a (t) = 9e^ B »V- B ««- 5 »- l S) 2s2H (20) 

the term — | (|) 2 5 2H being included to insure that E (a (t)) = 9. 

Eqs. (|TT|) and (fT%|) . or equivalently (fTTj) and (j20J), define our stochastic 
volatility model. 
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This stochastic volatility model will be used in the next section to derive 
an option pricing formula. However it is also useful to compute the modifi- 
cations implied by the stochastic volatility on the probability distribution of 
the price returns. From (JTSj) one concludes that logo* is a Gaussian process 
with mean (3 and covariance 



tJ, (s, u ) = ^ ||s - u + S\ 2H + \u - s + 5\ 2H -2\s- u\ 2H } (21) 

This Gaussian process has non-trivial correlation for H ^ |. At each fixed 
time logcr t is a Gaussian random variable with mean (5 and variance k 2 5 2H ~ 2 . 
Then, 

/ S 1 n 1 f (log(7-/3) 2 ] 

p s (a) = -ps (log a) = exp - (22) 



therefore 



with 



S T \ f°° ( S T \ 

log J = J d(J PS 0) Pa \log —J (23) 



,, s T \ i \ (io«(t)-(/"-*)(r-<)) : 

p a log — = == exp ' 



ft; ^trx 2 (T-t) I 2a 2 (T-t) 

(24) 

One sees that the effective probability distribution of the returns depends 
on the observation time scale 5. This is a pleasant feature of this stochastic 
volatility model, which contrasts, for example, with GARCH models which 
describe well the volatility at a given time resolution, but fail to account for 
the S dependence. 



3 Option pricing. Risk-neutrality approach 

Assuming risk neutrality the value V (St, <7t, t) of an option must be the 
present value of the expected terminal value discounted at the risk-free rate 

V(S t ,a t ,t) = e~ r ^ J V(S T ,cr T ,T)p{S T \S t ,a t )dS T (25) 
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V (St, cr T , T) = max [0, S — K] and the conditional probability for the ter- 
minal price depends on S t and a t . As in Hull and White [7j, we make use of 
the relation between conditional probabilities of related variables, namely 

p(S T \S t ,a t ) = J p (S T \ S t , log a) p (log a\ log a t ) d (log a) (26) 



logcx being the random variable 

1 ' T 



loga = f~t J iogasds ( 27 ) 
Then Eq. f|25|) becomes 

V(S t ,a t ,t)= fc^S t ,e^,t^p(hg^\loga t )d(\^a) (28) 

C^e^t) = J e^ T ^V(S T ,a T ,T)p(S T \S t ,hg^)dS T (29) 

C (s t , e log<J , tj being the Black-Scholes price for an option with average 
volatility e logcr , which is known to be 

C (S t , a, t) = S t (a + b)N (a, b) - Ke'^^ (a - b) N (a, -b) (30) 

with 

(31) 



» = i(S5 +,V!rri ) 

and 



1 Z" 00 2 

N{a,b) = -= dye- v -^ a+h)2 (32) 

v27r y_i 



To compute the conditional probability p (log a\ log 07) we recall that from 
our stochastic volatility model 

loi^ = log a t + j^— t J - 5 ds J (dB H (r) - dB H (r - 5)) (33) 

Notice that, because we want to compute the conditional probability of log a 
given log a t at time t, a t in Eq. (p£j|l is not a process but simply the value of 
the argument in the V (S t ,<J t ,t) function. 
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As a t— dependence process the double integral in is a centered Gaus- 
sian process. Therefore, given loga t at time t, loger is a Gaussian variable 
with conditional mean and variance 

E {logcr| logcr t } = logdi (34) 

with 

h = {2H + 1 ) {2H + 2) + 5 ) 2H+2 + {T-t- 5) 2H+2 — 2 (T — f) 2 ™ - 25*™} 

(36) 

J a = — — {2 (T - t) 2H+1 — (T — t+ 5) 2H+1 -(T-t- 5) 2H+1 ) 
2H +11 J 

Because in general 5 <C (T — t) one obtains 

« 2 - ^ (l - (2" - 1) ( j^T" I (37) 



Finally 



and from 



n 11 \ 1 J -(logo- -logo*)' 1 

P (log <t\ log a t ) = — ==- exp <^ — * — ^ (38) 



/oo 
dfC($,e*,t)p(f|log<r t ) (39) 
-00 



one obtains 



V (S t , a t , t) = S t [aM (a, a, b) + bM (a, b, a^-Ke'^'^ [aM (a, a, -b) - bM (a, -b, a)} 

(40) 

M{a,a,b) = / dy dxe~^ e'^^-) (41) 

2vra J_ x J 

1 [2 f 00 , c f ax 

ax rertc = — 



4a V 7T J ax + I \ V2 y/2. 
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T-t 0.5 S ,k 



T-t 0.5 g/ K 




T-t 0.5 S/K T-t 0.5 



Figure 2: Option price and equivalent implied volatility in the risk-neutral 
approach to the stochastic volatility model 

as the new option price formula (erfc being the complementary error function 
and a and b being defined in Eq. (|31j) ). with a replaced by a t . 

In Fig. 2 we plot the option value surface for V (S t , &t, t) in the range T — 
t e [5, 100] and S/K G [0.5, 1.5] as well as the difference (V {S h a t , t) - C (S t , a t , t)) /K 
for k = 1 and k = 2. The other parameters are fixed at o = 0.01, r = 
0.001,5 = 1,H = 0.8. The lower right panel shows the implied volatility 
surface corresponding to V (S t ,a t ,t) for k = 1. Notice that computationally 
it is perhaps more convenient to use Eq. (and an optimized Black-Scholes 
routine) than the more elegant Eq. ()40|) . 

4 An option pricing equation 

Whenever the actual drift of a financial time series can be replaced by the 
risk-free rate we are in a risk-neutral situation. However, in stochastic volatil- 
ity models this may not be a good assumption. In particular, an accurate 
option pricing formula should take into account the market price of volatil- 



11 



ity risk. Here we derive an option pricing equation which does not assume 
risk-neutrality. 

Because volatility is not a tradable security, a pure arbitrage argument 
cannot determine completely the fair price of an option. On the other hand, 
because of the fractional nature of the volatility process, volatility follows a 
stochastic process different from the one of the underlying security. There- 
fore, we cannot apply the reasoning that leads to uniform coefficients of 
the form (/ij — Aj<7i) in the first derivative terms of the option pricing equa- 
tion 1 . Hence, a first principles derivation, with clearly specified assumptions 
is required. 

As in Black-Scholes, we form a portfolio 

U(t)=V(S,a,t)-A(S,a,t)S t (42) 
From (|TT|) and (fTHJ) and choosing A (S, a, t) = |^ we obtain 

{f)V 1 r) 2 V 1 r)V k 

Consistent with our stochastic volatility model (JUJ), we have assumed the 
volatility to be uncorrelated with S and Eq . (J43|) follows from the application 

of the fractional Ito formula EH |2B|- Namely, if X t = Xf\ • • • , X t (n) ) 

with dxf } = a (t, to) dB { S (t) , then 

df (t, x t ) = %dt + j2 ^ dx P + E ^w Ci (t ' u) D * {Xt) (44) 

i i 

Df t (X t ) being the <\>— Malliavin derivative corresponding to the Xj: — process, 
defined by 

Dt )f X t (ui) = ton^X\ui + ej dsj <P(s,u)f(u)duj-X(iu i l^5) 
D* u (X t ) f (u) du 



Hi, Oi and Xi would be the drift, volatility and market price of risk for each process 
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(s, u) being the kernel 

<f>(s,u) = Hi(2Hi - 1) \s -u\ 2H ^ 2 (46) 

In (gnj) we are still left with the stochastic term er|^| (dB H (t) - dB H (t - 5)) 
and, because volatility is not a tradable security this term cannot be elimi- 
nated by a portfolio choice. Instead we may assume as reasonable to equate 
the deterministic term in dU (t) to 

k f)V\ 

rU (t) + v y±\ dt (47) 

where r is the risk-free return and, with v > 0, the second term is a measure 
of the market price of volatility risk (bigger risk, bigger return). We end up 
with 



dV n dV a 2 S 2 d 2 V k, 1TTr2H2 N dV rri2r2H , 2 d 2 V 

W + rS as + — as^ + I ( kHS - ") " a7 + Hks ° V = rV 

(48) 

as the general form of the option pricing equation consistent with the stochas- 
tic volatility model in Q. 

We now obtain an integral representation for the solution of this equation. 
With the change of variable 

S 

x = \og— (49) 
and passing to the two-dimensional Fourier transform 

V (t, x, a) = J J d<pdpF (0, p, a) e^ t+px) (50) 

we obtain 



Now, defining new constants 



(51) 



X (p) — 2Hk5' 2H - 2 

e( P A) = X 2 ( P )- r w^m (52) 

^2 ( n \ *P+P 2 
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and making the replacement 

F(a) = aXZt{(a) (53) 

Eq. (jHIl) reduces to a standard Bessel equation. Therefore the solution of 
is 

V (jt, x,a) = J j dpd^+^a^Z^) (C (p) a) (54) 

((a) being a Bessel function. The Bessel function will be a linear combi- 
nation 

Z i {Ca) = c l J^{C ! a) + c 2 N i {Ca) 

of a Bessel function of first kind and a Neumann function, with coefficients 
Ci and C2 to be fixed by the boundary conditions, which for call options is 

V (T, x, a) = max (x, 0) 

Eq.(JS2J) is an exact solution of the option pricing equation. How use- 
ful it might be will depend on the development of adequate computational 
algorithms for its numerical evaluation. 
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